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h— 5 1 Abstract 

CN We prove the universality of correlation functions of chiral complex matrix models in the 

microscopic limit (N — > do, z —* 0, iVz = fixed) which magnifies the crossover region around the 
origin of the eigenvalue distribution. The proof exploits the fact that the three-term difference 
0^ | equation for orthogonal polynomials reduces into a universal second-order differential (Bessel) 

■ equation in the microscopic limit. 
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The concept of universality in the theory of random matrices, or the independence of 
relevant quantities upon details of the potential function, is crucial in its application toward 
the level statistics of disordered physical systems. The universality of the macroscopic bulk 
two-point function in the large- N limit is discovered to hold as a direct consequence of 
the linear functional relationship (Cauchy inversion) between the potential V(M) and the 
large- N spectral density p(z) = limAr^oo (1/Ntv5(z — M)) 



Viz) r dw , . , . 1 r dw a 2 -z 2 V'(w) , . 
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0, and that p(z, w) ~ 5p(z)/SV(w). An alternative proof using orthogonal polynomials is 
found in ref.||. 

Another class of universality of a different origin, termed microscopic, has been anticipated 
whenever a quantity in concern is governed by the microscopic repulsion (anti-crossing) 
between energy levels dominantly enough to surpass the effect of a slowly varying potential. 
This type of (conjectured) universalities includes the appearance of the sine kernel for the 
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microscopic correlators in the bulk |4|,[| (proved in ref . ) , of the Airy kernel in the 
vicinity of a 'soft edge' (p(z) ~ y/a — z) J7|a (proved in ref.§), and of the Bessel kernel 
of a 'hard edge' (p(z) ~ 1/V° z ) [1>[[0 ^- The last and open problem has attracted 



considerable attention from condensed matter physics on spin-impurity scattering JTTJ] as 
well as from high energy physics on QCD chiral symmetry breaking and will be the 
central subject of this article. 

The problem can equivalently be formulated in terms of chiral complex matrix models as 
follows: Consider a matrix integral with a generic potential 



Z = f d 2N2 M exp {--tr V(M 2 ) ) , V(M 2 ) = ]T ^M 2k (2) 
J 12 J fe > x k 

where M stands for a 2N x 2N block hermitian matrix whose non-zero components are 
N x N complex matrices on the off-diagonals, 




M= . (3) 



A 'chiral' complex matrix model (or chiral unitary ensemble) is so called because of the 
invariance under the transformation 

C^UCV\ U,V eU(JV). (4) 

Since M anticommutes with 75 = ^ °^ ^ , each of its eigenvalues z always accompanies its 

mirror image —z in the spectrum. The repulsion between these pairs is expected to yield a 
region in the vicinity of the origin where eigenvalues avoid populating. Magnify this region 
by measuring the correlation functions in the unit of an average level spacing ~ 1/N, i. e. 
by substituting z = (/N. Are all the correlators independent of the potential in the limit 
iV —>■ 00, when ( is kept fixed? 

The answer has already been conjectured affirmative on several grounds [TSyTJ]. First, one 
naively expects that, in the vicinity of the origin, the potential can essentially be regarded 
as a constant. Therefore it would not affect the microscopic correlation due to the level 
repulsion, except via the average level spacing used as a unit, which is determined by global 



balancing. Ref.|T3[ noticed that one and the same Bessel kernel emerges from two ensem- 
bles with distinct potentials, Gaussian V = M 2 and an infinite wellB. This suggests that 



the Bessel kernel holds for a generic potential. Ref. [Tl has calculated the one-point func 



tion for V = M 2 + g M 4 up to O^g 1 ), which again supports the conjectured universality. 
We shall give a rigorous proof of this conjecture by exploiting the fact that the three-term 
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difference equations for orthogonal polynomials, characteristic of one-matrix models, uni- 
versally reduces the Bessel equation in the above mentioned microscopic limit. 



The partition function for a chiral complex matrix model (|2|) is expressible in terms of 
the component matrices as well as of the eigenvalues after integration over the angular 
coordinates (U, V) G XJ(N) x V(N)/XJ(l) N } 



d N2 C^ d N2 C e~ NtrV ( ctc ' ) 

00 TV 00 N 

= J n K 2 e- Nv ^) n k 2 - 4 = J n K e-^<*>) n ia, - a/ . ( 5 ) 
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The above expression can be interpreted as a positive definite hermitian matrix model in 
H = C^C whose eigenvalues are Ai, • • • , Ajv > 0. In other words, the problem reduces to 
finding a set of orthogonal polynomials P n (A) over the semi- infinite interval [0, 00), 

00 

J dX e~ NV ^ P n (A) P m (A) = h n 5 nm . (6) 


We normalise them such that P n (0) = 1 for later convenience, 

P n (X) = l + ---+ Pn X n . (7) 

Here it is assumed possible to choose this normalisation, which will prove equivalent to 
the ansatz that the origin be included in the support of the large- N spectral density of M. 
Then the recursion relation for these P n 's reads 



AP n (A) = -g n P n+1 (A) + s n P n (A) - gVi-17— — -Pn-i(A) [Qn ~ ''" 



iin-1 V Pn+l , 

— ^ t X nm P m ( A) . 



The sets of unknowns {h n }, {q n }, {s n } are iteratively determined by |16 



1 

1 = -/ dX—{e- NV ^ P n (X) P n (X)} = NV'CX) nn h n , (9) 





00 



(10) 
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and eq.(||) at A = , 

= -q n + s n - g n _ x — (11) 



We can immediately eliminate s n 's using (|TTD to get 



XP n (X) = -q n \p n+1 (X) - P n (X) - — |^ (P„(A) - P n -i(A)) 1 . (12) 
In the following we need to know the asymptotic behaviour of q n and h n for 

TL 

n, N — > oo while — = t is kept fixed. (13) 
Eqs.(|), (0) and (pD tell us that they should behave asi 

Tl I TX \ 

In = q(-tt) + higher orders in — , N h n = h(—) + higher orders in — . (14) 
N n N n 

Then the matrix A and its powers are approximated to be 

Tl 

(A'W = q{^) k £ (-)' ( , ^ n ) S nm+t (15) 
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so that eqs.@ and (|ToD read 



V «M M =55rw (17) 



fc \ 



By eliminating h(t) out of the above two, we obtain an algebraic equation for q(t), 

5 = 



Eqs.Qlq) and (|T8D imply a universal relationship among total derivatives, 

dt = 2q d (~) + i dg = 1J-q d . (19) 



t g n and /i ra converge to smooth functions when the eigenvalues are supported on a single interval. 



Next we expand the rhs of the recursion equation flT2|) in terms of 1/n in the limit (|13"D , 



A P(n, N, A) = |i + £ ^ j P(n, iV, A) + higher orders in i (20) 

where the argument N in P(n, N, A) = P n (X) is to indicate explicitly the dependency via 
the coefficient in front of the potential. It equivalently reads (subleading terms suppressed) 

telling us that that the arguments of P appear only in the combinations t = n/N and x = 
N 2 X in the limit (|13D. The rescaled eigenvalue coordinate x is to be fixed finite hereafter, 
and is regarded as a parameter in the ordinary differential equation in t. Performing the 



change of variable 1 1— > u{t) = yq(t)/h(t), using the relationship (|TU|) and neglecting higher 
order terms in 1/n, eq.(^) reduces to the Bessel equation of zeroth order, 

(-— u^- + Ax ) P(u,x) =0. (22) 
\u au au J 

The general solution to it is a linear combination of Bessel and Neumann functions 

P{u, x) = c{x) J (2u^) + c'(x) Y (2u^) . (23) 

The integration constants (functions in x) are completely fixed by the boundary condition 
at t = n/N = (it(0) = 0), 

P(0,x) = P (X) = l, (24) 

to be c(x) = 1, c'(x) = 0. By substituting t = 1 into (|23"D and its t-derivative, we establish 
the following lemma for the asymptotic behaviour of generic orthogonal polynomials over 
the semi-infinite range [0, oo) which are normalisable to P n (0) = 1: 



JffL P "(^ = J o{ 2u WVi), (25) 



i™ N { p "0 - Pw - i( ^0 = -iw> Ji (2u(1) ^ • (26) 

The parameters u(l) and q(l), through which the dependence upon the potential enters 
the asymptotic form of the orthogonal polynomial, have the following simple meaning. 
Comparison of eqs.(|16D and ( pT8|) at t = 1 with an explicit expression for the rhs of eq.([l|) 
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enables us to relate the parameters with the edge of the support of the large- N spectral 
density p(z) of the chiral complex matrix and its value at the origin as 

u(l) 



a = 2Jq(l), p( z = 0) = ^. (29) 
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It is easy to check that the critical condition p(0) = for the l-/2-cut transition WJft is 
equivalent to P/v(0) = for the conventional, monically normalised polynomials P n (X) = 
X n + • ■ -. Thus, under the assumption that the normalisation (|7|) is possible, the constants 
h(l) and q(l) are determined positive, due to the identification (p9|) valid for p(z) supported 
on a single interval. The positivity of the norm h is necessary for consistency, whereas q > 
signifies that the coefficients in P n (X) are alternating and eventually causes an oscillatory 
microscopic spectral density as it should. 

Now we recall the expression for the integration kernel Kn(X, p) associated with the eigen- 
value problem for the positive definite hermitian matrix H = C^C, 



K N (X,p) = e-T(vW+vM) Iv ^iS 

N £o h i 

= e _^(y(A)+y( At )) -QN-1 Pn(X)Pn-i(h) - -P/V-l(A)P/y(/i) 

Nh N ^ X-p ' { ' 



Here use is made of the Christoffel-Darboux identity. Plugging in the lemmata fl2"o] ) and 
(p6|), we obtain a universal form of the kernel (called the Bessel kernel) in the microscopic 
limit 



lim — Kat(— ^,-^t) = 

^ 70(2^(1)^)^(2^(1)^) - y/iJ! (2^(1)^.70(2^(1)^) 

ar-j/ 

The s-point correlation function (pn) of eigenvalues of H (M) is represented in terms 
of the kernel as 



a N (X 1 ,---,X s ) = (f[^-ti6(X a -H)) = det K N {X a ,X b ) (32) 
\ a =i N I ^ a ' b ^ s 

p N {z u --- ,z s ) = (t[ ^to8(z a - M)^ = \zi\---\z s \a N (zf r --,z^), (33) 
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respectively. Therefore all the formulae for their microscopic limits 



? (a*, • - • , x s ) = Jim ^(^, • • • , ||) (34) 

g(Ci, ■■■,Cs) = j^Piv(f , " " " , |) = ICil " " " IC.I ?(C?, • • • , C 2 ), (35) 

previously calculated for the Laguerre (in the if-picture) or chiral Gaussian (in the M- 
picture) unitary ensemble, hold universally. Namely, the spectral density of the chiral com- 
plex matrix model 



Pn(z) = (^ti6(z-M)) = \z\K N (z 2 ,z 2 ) (36) 
universally takes the form 

Q(0 = ICI (^o 2 (27rp(0)C) + ^(2^(0)0) (37) 

in the microscopic limit. It enjoys the matching condition between the micro- and macro- 
scopic (large- A/") spectral densities, 

Km g(C) = p(0). (38) 

(,^00 
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In this article we have exhibited a proof of universality of the correlation functions of 
chiral complex matrix models in the microscopic limit. Consequently all dependencies of 
correlators upon the potential appears only through a single and local (at z = 0) parameter 
p(0) as anticipated. The universality holds as long as p(0) > 0, i.e. the large-A^ spectral 
density is supported on a single interval. 



We have extended this strategy of taking the continuum limit of recursion equations, 
commonly used in the double-scaling calculations, to show microscopic universalities in a 
variety of matrix models which have been argued to be of relevance in QCD Details 
of these generalisations will appear elsewhere fi"9 
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